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1. INTRODUCTION 


To develop the next generation of high-speed civil transport (HSCT) aircraft, a significant 
reduction in pollutant emissions from the engine must first be demonstrated. Emissions that must be 
reduced include NOx emissions into the stratosphere during supersonic flight and the NOx, CO, and 
unburned hydrocarbons emissions near urban areas during idling and takeoff/landing operations. One 
candidate low-NOx combustor is a dual-stage combustor consisting of a fuel-rich primary stage, a 
quick quench zone, and a secondary fuel-lean burn stage. Significant NOx production appears to 
occur primarily in the quick quench zone during the unsteady mixing process. Numerical steady-state 
prediction models arc incapable of including the effects of unsteady mixing. This Phase I work 
demonstrates a novel unsteady mixing model that, when used in current numerical prediction codes to 
account for the mixing process, will provide the necessary correction, resulting in an improved 
analysis tool of great use to the engine designer. 

Extensive investigations have been carried out by various researchers (e.g., McVey and Kennedy, 
1979; Westmoreland et at., 1979; Bahr, 1980; Mularz, 1979; Cooper, 1979; Semerjian et a)., 1979; Tang 
and Churchill, 1981; Tacina, 1990) to identify and evaluate various advanced combustor designs that 
could result in reduced pollutant emissions. It is well-known that NOx production is kinctically 
controlled (an exponential function of adiabatic combustor temperature) and occurs due to the 
availability of free oxygen atoms in hot regions. Therefore, a low-NOx combustor design must 
minimize or eliminate free oxygen atoms in the hot regions without decreasing the combustion 
efficiency of the combustor and without increasing the emission of CO and unburned hydrocarbons. 
Various low-NOx combustor concepts — i.e., Lean-Premixed-Prevaporized (LPP), Rich Burn/Quick 
Quench/Lean Burn (RQL), and Direct Injection (Dl)-are being studied and their NOx emission 
characteristics have been determined over a wide range of inlet temperatures, fuel-air ratios, and 
pressures (Tacina, 1990). Current data indicate that the LPP concept, which involves premixing fuel 
vapor with air and burning at low temperature, has the lowest NOx emission. However, this concept 
has a major operational disadvantage due to its narrow stability limits. The DI concept, which 
involves direct injection of the fuel into the combustor, has shown promise in that it has increased 
stability limits and also produces low NOx emissions comparable to that produced by the LPP 
combustor. However, results for the DI concept so far have been obtained with gaseous fuels, and it 
is not clear whether the same characteristics will be observed with liquid fuels (Tacina, 1990). 

For aircraft such as the HSCT, which will have a wide range of operating conditions, the RQL 
combustor configuration has good potential for significantly reducing NOx emissions. An idealized 
RQL configuration is shown in Figure 1. In such a configuration, the primary combustion occurs in a 
fuel-rich environment. Downstream, air is injected to rapidly dilute the fuel-rich mixture exiting from 
the primary zone so that a uniformly mixed fuel-lean mixture burns in the secondary combustion zone. 
This concept, although more complex than the LPP and the DI concepts, attempts to minimize (if not 
completely inhibit) NOx formation in the hot primary stage by limiting or eliminating free oxygen 
atoms and in the fuel-lean secondary stage by reducing the combustion temperature. NOx will form 
in regions where the secondary air mixes with the hot fuel-rich mixture exiting from the primary zone 
(in the quench region), especially in regions where the fuel-air ratio becomes locally close to 
stoichiometric. In fact, current data appear to indicate that most of the NOx may be forming in the 
mixing region. However, the exact mechanism of NOx production in the highly unsteady mixing zone 
has not yet been clarified. The RQL concept has good stability characteristics due to its fuel-rich 
zone and can be used with current jet fuels as well as alternative fuels containing large amounts of 
fuel-bound nitrogen (Tacina, 1990). Thus, if an effective technique to rapidly mix the secondary air 



with the fuel-rich mixture can be determined, this concept may become practically feasible. 

At present, there is an extensive experimental program underway to evaluate Iow-NOx combustors 
such as the LPP and the RQL (Ott, 1990). Detailed measurements inside and immediately downstream 
of the quench zone will be required to determine the efficiency of the rapid mixing concepts. Experi- 
mentally studying various types of rapid mixing concepts such as swirl jets and jets in crossflow can 
be expensive, and unless the characteristic length and time scales involved in these unsteady mixing 
concepts are determined and their relevance to NOx production quantified, the experiments alone 
may not provide the necessary guidelines to build production engines. Typically, numerical predic- 
tion methods are used to complement experimental studies, both to understand the dynamics and to 
extend the experimental results to unmeasured or unmeasurable operating conditions. However, most 
numerical schemes currently used for prediction purposes are based on steady-state schemes (e.g., 
Correa, 1984; Correa and Shyy, 1987; Magnussen and Hjertager, 1977; Swithenbank et al., 1980; 
Nguyen and Bittker, 1989) and are unable to account for the highly unsteady phenomena occurring 
during rapid mixing. Some earlier attempts to include the effects of unsteady mixing (e.g.. Pope and 
Correa, 1986; Pratt, 1980; Butler and Pratt, 1986; Chen and Dibble, 1990) have resulted in somewhat 
improved predictions. In principle, simulation techniques such as large-eddy simulations (LES) could 
be used to understand and predict more accurately the unsteady mixing processes, but they arc 
impractical for engineering design analysis due to the substantial computational resources required. 
Therefore, there is a critical need to develop a more accurate yet simple model for unsteady mixing 
processes that can be used in conjunction with conventional steady-state methods so that the com- 
puted results are corrected for the unsteady processes occurring in the mixing region. The develop- 
ment and demonstration of such an unsteady mixing model are described in this report. 

A numerical approach that models unsteady flow fields must reflect two practical objectives: (1) 
the physics of the flow field must be modeled as accurately as possible to ensure reliable predictions, 
and (2) the numerical technique must provide the design engineer with information as fast as possible 
to improve and optimize the design. In reality, these two objectives are often in conflict, and thus a 
compromise is necessary. 

The mixing between the fuel-rich mixture and the secondary air is an unsteady process occurring 
in a highly turbulent environment where significant fluctuations in both the velocity and temperature 
field will occur. Velocity fluctuations will affect the turbulent transport and mixing of the various 
species, and the temperature fluctuations will affect the chemical kinetics rates leading to non- 
equilibrium production/destruction of species. In fact, the prediction of thermal NO production 
based on the assumption that the O atoms are in equilibrium with O a (Peters and Donnerhack, 1981) 
has been shown to be in error due to significant non-equilibrium concentration of O and other radi- 
cals in turbulent nonpremixed flames (Drake et al., 1987). 

Another physical consequence of secondary air injection into the fuel-rich mixture is that in 
regions where these two streams come into contact, a shear layer will form. Such shear layer regions 
are well-known to be susceptible to instability and will roll up into vortices. These vortices then 
undergo pairing and merging processes, resulting in large coherent vortices that subsequently break 
down into small-scale turbulence. Such structures are known to play a major role in entrainment and 
mixing processes, and their effect on mixing needs to be taken into account. For example, there could 
be regions where the local fuel-air ratio is near stoichiometry, and this could result in large NOx pro- 
duction if the temperature is high. Rapid mixing concepts typically attempt to break down these 
structures as quickly as possible so that turbulent mixing will be enhanced. However, there may still 
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be regions in the mixing zone where these vortices persist for some time, and so the characteristic 
length (or size) and time scales associated with their formation and breakdown need to be included in 
the unsteady mixing model. 

Studies of molecular mixing have shown that the mixing process generally consists of two steps: 
(1) the entrainment of the fluid from the two streams into the mixing region, and (2) the mixing of the 
fluids at the molecular level as they come into contact. Classical modeling of turbulent mixing in high 
Reynolds number flows typically assumes that turbulent convection dominates over molecular diffu- 
sion processes, i.c., the first step is the crucial one. However, recent studies at CalTech (for a sum- 
mary, see Dimotakis, 1989) have revealed that, for flows at the same Reynolds number, both the mix- 
ing rates and the distribution of the mixed fluids were quite different in fluids of widely different 
molecular diffusivity. This clearly demonstrated that, besides turbulent convection, molecular diffu- 
sion effects are also important at the small scales. In terms of nondimensional parameters, this implies 
that the effects of both the Reynolds number and the Schmidt number on the mixing rate need to be 
taken into account. When chemical reactions are also occurring, the Damkohler number effect will 
also need to be included. 

Kerstein (1988, 1989, 1990, 1991a, 1991b) took into account these earlier results and developed a 
novel model called the linear-eddy model, which treats separately the effects of both turbulent convec- 
tion and molecular mixing at the small scales. It was demonstrated that this model is capable of accu- 
rately describing many features of turbulent mixing and chemical reactions in turbulent shear layers, 
nonpremixed flames, and turbulent jets. This successful demonstration of the linear-eddy model has 
laid the basic groundwork to further develop the mixing model for application to practical problems 
such as the prediction of NO production in the quench zone of the RQL combustor. In this Phase I 
study, a configuration-invariant mixing model has been developed so that it is applicable to a wide 
variety of flows. Furthermore, this mixing algorithm has been developed so that it can be imple- 
mented within the framework of well-used prediction codes to account for the unsteady mixing 
effects. This coupling issue will be addressed in Phase II. 
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2. PHASE I TECHNICAL OBJECTIVES 


The overall Phase I technical objective was to demonstrate the potential of the linear-eddy 
unsteady mixing model to take into account the characteristic mixing lengths and time scales within a 
mixing zone. To demonstrate its capability, this model was used to study a mixing situation similar to 
that which would occur in the quench zone of an RQL combustor. 

The specific technical objectives of the Phase 1 study were: 

• To formulate the linear-eddy model to characterize the unsteady mixing between the fuel-rich 
primary mixture and the secondary air stream: 

• To implement a reaction mechanism within the linear-eddy mixing model; and 

• To demonstrate the ability of the linear-eddy model to predict NO production in the mixing 
zone. 

As described in this report, all the above technical objectives were successfully achieved. The results 
of this research will be presented in the near future (Mcnon ct al., 1992). 
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3. FORMULATION OF THE MIXING MODEL 


The basic model and its application to a wide range of turbulent flows with and without chemical 
reactions has been described by Kerstein in a series of papers (Kerstein, 1988, 1989, 1990, 1991a, 
1991b). The motivation for these studies was to develop a configuration-invariant approach to study 
turbulent flows. Kerstein’s results showed that a set of laws or models governing the processes of 
entrainment, turbulent mixing, and thermochemistry can be developed to provide a self-contained pic- 
ture of the overall fluid mechanical and thermochemical structure of turbulent reacting flows. His stu- 
dies further showed that the same model is capable of predicting, for example, the observed differ- 
ences between shear layer mixing and thermal mixing layer in grid turbulence (Kerstein, 1988, 1989). 
In this Phase I study, this model has been extended to study unsteady mixing and thermochemiral 
processes similar to those occurring in a mixing zone. Multispecies mixing, propane-air nonpremixed 
flame and hydrogen-air nonpremixed flame problems were computed to demonstrate the capabilities 
of the mixing model developed in this research. 

3.1 The Linear-Eddy Unsteady Mixing Model 

The basic idea of the original formulation is to treat separately the two different mechanisms act- 
ing to describe the evolution of a scalar(s) (chemical species) in a specified domain. These two 
mechanisms are molecular diffusion and turbulent convective transport. For flows in which finite-rate 
kinetics is occurring, some additional mechanisms need to be incorporated, as will be discussed in the 
Section 3.2. 


3.1. 1 Molecular Diffusion 


The first mechanism is molecular diffusion, which is implemented by the numerical time integration 
of the diffusion equation 


dTk 

dt 


+ Dk 



(I) 


where Y k is the &th species mass fraction, W k is the chemical produciion/destruction term for the ktn 
species, and D k is the fcth species diffusion coefficient. In principle, this diffusion equation could be 
solved in the full three-dimensional form. However, the strategy here is to resolve all the relevant 
length scales involved in the local diffusion and transport processes. Note that it is in the small scales 
where mixing and chemical reactions occur; thus, resolution of the small scales is important. Conven- 
tional time-accurate simulation techniques (e.g., direct numerical simulations) attempt to resolve all 
the small scales numerically, resulting in prohibitive computational cost. Also, the resolution possible 
using current supercomputers only allows simulations of flows at low Reynolds numbers, which are 
not indicative of realistic systems. Even when subgrid models are used in LES techniques to increase 
the Reynolds number, the computational cost is still quite high. 

Therefore, to keep the computational cost reasonable and to resolve all the small scales taking part 
in the unsteady mixing, the domain is restricted to one dimension; hence the term linear-eddy. Ker- 
stein (1988, 1989, 1990, 1991a) applied this model to various types of turbulent shear flows and 
obtained good agreement with experimental data. Note that although Equation (1) is solved in one 
dimension, the information on the volume of the cells is built into the algorithm implicitly. This 
aspect will be described later. 
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For the purpose of discussion, consider a reaction mechanism between two species, A and B, such 
that 


A + B -* P (2) 

More complex reaction mechanisms typical of propane-air and hydrogen-air nonpremixed flame 
combustion have been studied in Phase I, as discussed later. The reaction rate for the mechanism in 
Equation (2) is 

W = ki[A\B] ( 3 ) 


Thus, for this reaction mechanism, the Fickian diffusion equation (1) becomes in one dimension 


3Y a 

dt 


- D 


ab 


d 2 Y a 

dx 2 


di\ 

dt 



-kY a Y b 


( 4 ) 


where k is the normalized rate constant. 

The application of the linear-eddy model to a particular flow configuration depends both on the 
configuration and on the issues that need to be addressed. In the earlier formulations (Kerstein, 1988, 
1989) the linear-eddy dimension (the computational domain) was chosen as a transverse (y-direction) 
line co-moving with the mean streamwise flow, and the streamwise development of the mixing region 
(x-direction) was obtained by relating it to the temporal evolution in the diffusion equation by the 
relation x = l/t, where U is the convective velocity. In a more recent study of turbulent nonpremixed 
flames (Kerstein, 1991a), it was shown that the choice of the streamwise direction as the linear-eddy 
computational domain is more advantageous, especially for the study of the axial structure of the flow 
based on the fluxes of mass, momentum, etc., through transverse planes. 

Consider the configuration shown in Figure 2, which shows a fuel jet exiting into a mixing zone 
This is a test geometry that corresponds to the mixing zone in the RQL combustor and other 
geometries similar to diffusion mixing and combustion in many experiments. Thus, demonstration of 
the mixing model in this geometry should provide the necessary validation for this Phase I study. 

If the x-direction is chosen as the linear-eddy domain, and if that domain is further subdivided 
into small cells, then each cell would represent a control volume, which is determined by the cell 
width (Ax), and a radial zone extending from the centerline to the nominal jet radius (r y ). The 
volume of each cell of width Ax is then given by V(x) = irr y 2 Ax. For the case shown in Figure 2, the 
jet radius is given by the relation 


r y (x) = C,(x - x 0 ) (5) 

where x 0 is the virtual origin and C, is an empirical constant which ranges from 0.07 in the forced- 
convection limit to 0.10 in the natural-convection limit The choice of Ax is dependent on the smal- 
lest characteristic length scale (or eddy size) that needs to be resolved. Kerstein (1991a) chose 
Ax = L(x)/n, where L(x) is the integral length scale, which is given by L(x) « 2r y (x), and n is the 
resolution factor, which is constant for a given computed realization. Typically, n = 50 was chosen 
for this study based on the analysis by Kerstein who found that this resolution is sufficient for many 
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of the earlier studies of turbulent mixing. The importance of the resolution will be addressed in more 
detail in the next phase. Using the definition of L(x ) given above, the volume of each cell becomes 

V(x) = -f L(x)*Ax . 

4 

Since we are interested in the downstream evolution of the mixing region, the spatial evolution of 
the jet needs to be incorporated. From experimental data, the spreading of the jet can be given as 

= 1 + cx (6) 

where d 0 is the momentum diameter of the jet, which is defined as d 0 = ^oVpo/pT and c is a constant 
denoting the spatial spreading rate; typically, c = 0.14 for free jet. Here, x = x */d& where x* is the 
dimensional streamwise distance. Also, p 0 and p a are respectively, the density of the primary fuel 
stream at the inlet (i.e., at x - 0) and the ambient (air) density. Note that, with this definition of 
L(x), the cell width Ax and the volume V(x) increases in the downstream distance. The increase in 
the volume reflects the entrainment of air into the mixing region. For the RQL combustor, the free 
jet spreading rate may not be applicable since mixing will be occuring in a confined region. There- 
fore, simulations with different values of the spreading rate parameter, c were performed to address 
the effect of modifying the effect of air entrainment. Some of these results are described in this 
report. 

In addition to the specification of the spreading rate of the mixing layer, it is also necessary to 
decide the total streamwise extent of the mixing layer that must be simulated. Typically, this distance 
is specified as the distance at which self-similarity is achieved. Experiments suggests that a region 
x /d 0 « 200 needs to be simulated since the premixed flame length typically can be in the range of 
x/d 0 = 100-140. In all the present calculations, the spatial extent of the mixing layer was divided 
into 25 blocks, and each block was then divided into 50 cells. Thus, the total number of cells used for 
the computation was 1250. Note that, since the integral length scale increases with downstream dis- 
tance, the size of the blocks increases in the streamwise direction. With the resolution used in this 
study, and with the spreading rate given by equation (6) (with c = 0.14), the total length of the mixing 
layer exceeds x/d 0 = 200. This was sufficient for most of the nonpremixed flame combustion prob- 
lems studied here. The effect of increasing the resolution and the streamwise extent of the mixing 
region have not yet been investigated, since this would increase the computational cost of each simu- 
lation. However, the calculations discussed in this report clearly demonstrate the capability of the 
numerical scheme, and the resolution used is considered sufficient for the Phase I validation study. 

Once the cell size (Ax), the spreading rate, the total number of cells, and the initial conditions (to 
be discussed later) are specified, the diffusion mechanism [Equation (1)] can be solved in a purely 
deterministic fashion using standard finite-difference schemes. 

5.7.2 Turbulent Convection 

The second mechanism in the linear-eddy model is the turbulent convection process. This is incor- 
porated in the model by a random process that rearranges the fluid element along the linear direction 
and is carried out (subject to some constraints) while the diffusion process is going on. This process 
causes a random walk of the fluid elements and introduces a discontinuous (turbulent) fluid motion, 
representing small-scale convective stirring. Each event involves permutation of the cells of the spa- 
tially discretized concentration field and is confined to a finite segment of the total computational 
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domain. This finite domain may be viewed as the size of the eddy that is causing the turbulent stir- 
ring. To bring the physics of turbulent transport into this process, the frequency of these events 
corresponding to a given range of eddy sizes is determined such that the principal scaling laws 
governing turbulent transport are satisfied. 

To implement this process of turbulent transport, two different rearrangement processes have been 
developed by Kerstein (1989, 1991a). These are the block inversion and the triplet map. The effect of 
both these processes are graphically described in Figure 3. The triplet map (Figure 3b) is used in all 
the calculations discussed in this report, and is described below. 

Figure 3b shows the effect of the triplet map in the continuum limit. Essentially, the size of the 
eddy (as described below) is first chosen from a probability distribution of eddy sizes. For this dis- 
cussion, let us choose an eddy size of length /. Then the scalar field within the chosen segment / is 
compressed by q factor of three, thus tripling the scalar gradient within the segment. The original 
scalar field within the segment / is then replaced by three copies of this compressed field, with the 
middle copy reversed. The resulting scalar field is shown in Figure 3b. This discrete mapping is 
defined so as to recover this rule in the continuum limit while satisfying species conservation exactly 
in the discrete implementation. Since in the present case the cell volumes are increasing in the down- 
stream direction, some care is required to take into account the unequal volumes of the cells. 

To describe the actual implementation of the triplet map, let us assume that the segment chosen of 
length / is in the discrete formulation made up of 12 cells. Let / consist of cells that are numbered 
1,2,3,4,5,6,7,8,9,10,1 1,12. When the triplet map is used to compress and rearrange the cells, the result- 
ing new sequence will be 1,4,7,10,11,8,5,2,3,6,9,12. The cells themselves are not physically moved; 
rather, the species in the cells are exchanged in this manner. 

The triplet map can be understood on both an intuitive and a theoretical basis. Intuitively, consider 
the effect of a single clockwise eddy on the scalar field that initially has a uniform gradient in the 
field. This is shown in Figure 4a. The effect of the eddy will be to distort the scalar gradient in the 
field, which would then take a form similar to that shown in Figure 4b. Comparing this figure with 
the effect of the triplet map (Figure 3b), it can be seen that the mapping procedure essentially reflects 
this phenomenon. More details of the theoretical considerations behind the triplet map are given in 
Kerstein (1990c). 

As noted earlier, to carry out the triplet map, the size of the eddy / must be first determined. The 
sizes of eddies present in the flow was assumed to range from the Kolmogorov microscale r? to the 
integral length scale L (i.e., r)<l <L) and the eddy size / was chosen randomly from a power-law dis- 
tribution /(/) within this chosen range. 

Once the range of eddy sizes is known, a single inversion event involves first choosing a eddy size 
/ from the pdf /(/) (which needs to be determined) and carrying out the spatial inversion described 
above. Note that /(/) varies with both x and /. A sequence of such inversions must be carried out so 
that all allowable block sizes (or eddies) within the chosen range take part in the process. The rate at 
which this inversion takes place is determined based on the local values and therefore varies both 
with time and with the spatial location in the mixing region. This inversion rate is specified by a rate 
parameter A, which has the unit of (length x time) -1 . This inversion rate is determined by first deter- 
mining the turbulent diffusivity of the random walk process for each eddy. The turbulent diffusivity 
of a random marker is then the sum of the contributions due to all the eddies taking part in the inver- 
sion events. Without going into detail (see Kerstein, 1989), the total diffusivity associated with the 
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random walk of all eddies up to a size l 0 (where r\ < / 0 < L ) is 


Did) = flrlVm 

n 6 


( 7 ) 


Using Kolmogorov scaling, Re^^do/ri) 3 ^ 4 , Kerstein derived an expression for /(/), which can be 
written as 


/(/) 


5 1 

3 L (L/r ,)*' 3 - 1 



( 8 ) 


Using this expression in Equation (7), the value of A can be determined by setting l 0 = L and 
D, = i/Re t . This is given as (Kerstein, 1989) 


A _ 54 t'R e L r(EA?) 5/3 ~ 11 
5 L 3 1 - ( r)/L)*> 3 


(9) 


For high Reynolds number flows, L » rj, the leading-order approximation to Equation (9) is 


A = 


54 


15/3 


( 10 ) 


In the algorithm, the full expression [Equation (9)] is used. The time scale associated with the ran- 
dom stirring events, r„ is related to the rate parameter A by the simple relation r, = 1 /XL. 

This completes the description of the basic model. To summarize, the unsteady mixing in a speci- 
fied domain is modeled by two processes that work together but which are modeled separately. 
Molecular diffusion updates the concentration fields deterministically based on Fickian diffusion 
within a linear domain in which all the relevant small scales are resolved. The typical time scale for 
this diffusion process is related to the cell spacing, Ax , and the diffusion coefficient D as 
id = Ax 2 /D. Chemical heat release within each cell in the linear domain can also be included in this 
diffusion process by solving the diffusion equations for k species with the full chemistry model for IV k . 
However, the finite-rate kinetics of interest in this study is a multispecies mechanism (described in the 
next section) with large production/destruction terms for W k (and thus, small chemical times, r ch ). 
This can make the diffusion equation numerically stiff. Therefore, we will implement the finite-rate 
kinetics in a different way within the linear-eddy, as described in the next section. 

While the molecular diffusion process is going on, it is punctuated by randomly occurring rear- 
rangement events, representing the turbulent convective stirring process. Thus, stirring occurs at a 
frequency A, as determined by the physical scaling laws described above, and results in the rearrange- 
ment of the species within each eddy of size /. Typically, the time (frequency) at which the inversion 
events take place will be different than the time-step required for integrating the diffusion equation. 
In most cases, r, « r d , and thus a series of inversion events have to be carried out with a characteris- 
tic time (r,) between each diffusion step. For example, if r d « nr„ then n inversion steps have to be 
carried out before the next diffusion step can take place. This is taken into account in the algorithm. 
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3.2 Application-Specific Processes 

The basic processes, molecular diffusion and turbulent mixing must be further supplemented by 
application-specific processes that characterize the type of unsteady mixing being studied. By 
parametrically varying the important parameters of the application-specific process, the effect of dif- 
ferent unsteady mixing strategies can be studied. In Phase II, these application-specific processes will 
be handled differently since the mixing algorithm will be incorporated with steady-state prediction 
codes. However, the physics involved in the application-specific processes will not be much different 
Some of the issues are described briefly below. 

Figure 2 shows a characteristic case in which a fuel jet and an air jet enter a mixing zone. A vari- 
ant of this type of mixing process is the more complex shear flow due to jets in crossflow. Such com- 
plex mixing scenarios will be considered in the next phase. For the configuration shown in Figure 2, 
there are three application-specific processes that need to be modeled. They are fuel inflow, air 
entrainment, and streamwise motion of the mixture. Each of these processes is described below. 

3.2.1 Fuel Inflow 

For the linear-eddy application, the streamwise direction is divided into cells of width Ax and 
volume F(x), as described above. Fuel enters the first cell of volume K(0) at a constant mass flow rate 
m 0 from the upstream state. The state conditions -e.g., the temperature, pressure, density, and the 
species mass fractions in the fuel mixture at the entrance to the mixing zone -need to be specified. In 
all the calculations described here, this state is specified. For the RQL combustor, this state may be 
the exit conditions from the primary zone. 

3.2.2 Air Entrainment 

The process of air entrainment is a major facet of this model. The entrainment process is deter- 
mined based on the length and time scales associated with the corresponding physical process. Thus, 
by proper modification of the entrainment law, the effect of different air injection strategies could be 
studied. Although this issue will be addressed in more detail in the next phase, some preliminary 
study of modifying the entrainment process was carried out in Phase I, and is described later. 

From previous experimental observations, it is clear that the air will be entrained throughout the 
flow domain; however, it is not entrained in a continuous manner. Rather, it appears that the air is 
entrained in finite parcels, and that and the size of the parcel directly influences the mixing process. 
From earlier experiments (e.g., Dahm and Dimotakis, 1987; Mungal and Hollingsworth, 1989; Mungal 
and O’Neil, 1989) it appears that the size of the parcel is comparable to the jet diameter. This infor- 
mation is included in the model by setting the volume of the entrained air parcel to be 
V E (x) = (x/4 )L\x). Note that, since L(x) is an increasing function with x, the volume of air 
entrained will also increase with streamwise distance. 

Thus, the entrainment of air parcels of volume V E (x) will occur throughout the domain. The fre- 
quency of the entrainment events in the interval (x, x + dx) is given by A E (x)dx, and is determined 
from the local entrainment rate from the relation 

= A £ (x)^K e (x) (11) 
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where m* is the mass flux of the air crossing any x plane. Based on the exit conditions at the fuel jet, 
the mass flux from the fuel jet m,{ 0, t) = m /0 (0 will be known. To obtain the mass flux of the 
entrained air at any x plane, m a (x), the use of entrainment laws is necessary. For example, the sim- 
plest law is given by Ricou and Spalding (1961) which was modified by Becker and Yamazaki (1978) 
to obtain 


dm* 

dx 



( 12 ) 


where p = m/V, and V and C are, respectively, the volume and momentum flux crossing any x plane. 
For this generalized law, C E is an empirical coefficient that ranges from 0.32 in the forced-convection 
limit to 1.84 in the natural-convection limit. 

The momentum flux G(jc) is determined from the local momentum balance according to the rela- 
tion 


% = *C 9 g( Poo -p)b 2 (13) 

where g is the gravitational acceleration, C g =\ 1.8 is an empirical coefficient, and b is the nominal jet 
radius defined earlier in Equation (5). 

For specific modeling of different mixing concepts, the entrainment laws need to be determined 
either empirically or from theoretical/experimental studies. In all the studies described in this report, 
Equation (13) is integrated with the inflow condition, G(0) = C 0 , where, G 0 = (x/4)p 0 Uo 2 ^o 2 where, 
subscript 0 denotes the incoming state of the primary fuel stream at the inlet. This result is used in 
Equation (1 1) to obtain the entrainment rate A E (x). 

The entrainment events occurring at the frequency A E are taken to be statistically independent 
(similar to the inversion events), and the implementation process is chosen to be similar to the tur- 
bulent stirring process described earlier. Thus, the epochs and locations of the entrainment events are 
selected by an algorithm that parallels the stirring algorithm described earlier. 

Note that when the entrainment events are taken to be statistically independent, the effect of 
coherent vortical motion cannot be included in the model. For accurate prediction of combustion in 
shear flows, the effect of coherent motion must be included since it is known to exist in shear flows 
and its effects on mixing and entrainment processes are well-documented. As noted before, however, 
organized entrainment patterns are configuration-specific; thus, to include their effect on the overall 
entrainment process would require specific input depending on the type of shear flow being studied. 
The exact modifications to the entrainment law necessary to include the effect of coherent vortical 
motions in the mixing layer will be considered in more detail in the next phase. 

3.2.3 Streamwise Motion 

When fuel enters the first cell, it must displace fluid in subsequent cells to maintain constancy of 
the cell volume. This displacement induces a streamwise flow. In some respects, this reactant-feed 
mechanism resembles earlier plug-flow studies (e.g., the coalescence-dispersion model of Pratt, 1980). 
However, there are two major differences. First, the radial spreading of the jet is included in the 
present model by virtue of the streamwise increase of the cell volumes (see Figure 2). Second, 
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although the fuel will enter at only one location (the exit from the primary combustor), the air will be 
entrained throughout the strcamwisc direction, as shown in Figure 2. Both these realistic mechanisms 
are included in the present model, whereas they arc missing in the plug-flow mixing models. 

3.3 Implementation of Reaction Mechanism Within the Linear-Eddy Model 


As noted in the above section, conceptually, there is no problem to include the finite-rate kinetics 
within each cell of the linear-eddy during the molecular diffusion process. However, the small time 
scales associated with the finite-rate kinetics will make the time-step very small and increase the 
overall computational cost. To reduce the computational cost, a more practical approach is used. Fig- 
ure 5 shows how the finite-rate kinetics algorithm fits into the linear-eddy algorithm. 

In this study, three different mixing and reaction mechanisms were implemented to demonstrate 
the capability of the mixing model developed in this research. Here, we will describe each of the 
approaches. 

3.3.1 Multispecies Mixing 

For a series of calculations, the effect of multispecies mixing processes in the jet was studied. For 
this demonstration, the primary stream or the fuel jet was chosen to consist of five species. Most of 
the calculations were performed using species 0 2 , N 2 , O, N, and NO with specified mass fractions at 
the inlet (x = 0) and with a specified inlet temperature. Some simulations were also carried out for a 
mixture that mimiced the distribution exiting from the RQL combustor. The air stream was assumed 
to contain only two species, 0 2 and N 2 , again with specified mass fractions and temperature. The 
mixing and diffusion process between the species from the fuel jet and the entrained air stream was 
then computed for various fuel temperatures, fuel jet velocities, and initial species distributions. The 
diffusion process was modeled by solving Equation (1) for Y k where k = 1,5. In some simulations, the 
effect of modifying the entrainment process was also studied. No reaction was assumed to occur for 
most of the mixing studies (IF* = 0). However, since the initial distribution contained both O atoms 
and N 2 , it was possible to estimate the formation of NO using a global NO production mechanism 
given by Chen and Kollman (1991). This mechanism can be used to estimate the formation of NO (in 
gm/cm 3 - s) by 


W NO = 2k J lN 2 ){0]M NO (14) 

where A/, is the ith species molecular weight, [C] denotes the concentration of species C (moles/cm 3 ), 
and kf = 1.84x 10 u exp(-38370/T) cm 3 /mole - s. It must be pointed out that this NO production 
rate was determined by Chen and Kollmann (1991) for the H 2 -air combustion mechanism (described 
below) and is used here only to demonstrate the effect of mixing on the NO production mechanism. 
The effect of finite-rate kinetics and the associated NO production was studied in more detail using 
the H 2 -air combustion mechanism as discussed below. 
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3.3.2 Propane- Air Nonpremixed Flame 

For this study, the earlier study by Kerstcin (1991b) was repeated using the new code. To study 
this type of flame, a reduced mechanism can be used in terms of the mixture fraction $, which at any 
x location is given by 


«*> = 


SSSL 

m(x) 


(15) 


where m(0) is the mass flux at the entrance (the fuel jet). Thus, the inflow condition (at x = 0) for 
this case is given by £ = 1 . 

The thermochemistry of the propane-air flame is treated by assuming equilibrium chemistry and 
adopting the Shvab-Zeldovich assumptions so that the local density depends only on the local mixture 
fraction, p = p(0 . Thus, the axial variation of the density can then be determined once $(x) is deter- 
mined. The relation between p and $ for propane-air combustion based on equilibrium chemistry 
(Gordon and McBride, 1971) is determined in the range 0 < $<0.15. For larger values of $, soot 
begins to form and thus invalidates the equilibrium assumption. The functional form is 


ASL = 

Poo 




J__Li 

k’l-G 


-1 


(16) 


where $* denotes the stoichiometric value (for the propane-air nonpremixed flame, $* = 0.0601), and 
T f is the adiabatic flame temperature (here, T f = 2283 K). This expression is valid for $ < 0.15; for 
0.15 < $ < 1.0, the density is obtained by extrapolation. 

To determine $(x), the diffusion equation (1) is solved with Y k replaced by $. Thus, the diffusion 
equation becomes 


21 m 

dt 


D 


£L 

' dx 2 


(17) 


The molecular transport coefficient D ( in Equation (17) and the molecular viscosity v are updated in 
each cell at each time based on a T 1 ’ 7 dependence, i.e., D(T ) = CT 1,7 and v = ScD, where. Sc is the 
Schmidt number which is a specified input (typically. Sc = 0.406). To obtain the temperature varia- 
tion, we used the equilibrium relation for density to get 


no. = i + 


min 

[j_ ±12-} 

T 

1 oo 

T 

* oo 




(18) 


Finally, an estimate of NO production in the propane-air flame by the thermal mechanism (the so- 
called Zelodovich mechanism) can be estimated. Here, we use the approximate expression of Peters 
and Donnerhack (1981) for the thermal NO production rate Sso per unit mass. This relation is 
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1/2 


(19) 


Sno 


BM ho 


I*L 

Ms t 


pro, 


exp(-<?/D 


where B = 5.74* 10“ (cm 3 /mole) 1/3 /sec and Q = 66,900. This expression was obtained by assuming 
that O atoms are in chemical equilibrium with 0 3 molecules. 

This expression has been criticized by Drake et al. (1987) because significant superequilibrium 
concentrations of O and other radicals occur in turbulent nonpremixed flames. They estimated that 
the prediction by Equation (19) is off by a factor of 2.5. 

To determine the NO production rate S NO for the propane flame, the terms in Equation (19) arc 
rewritten in terms of £. To obtain the mass production of NO, the NO rate in each cell is summed 
over all linear-eddy cells and divided by m(0) to obtain the emission index. 


3.3.3 Hydrogen-Air Nonpremixed Flame 

For the propane-air nonpremixed flame problem described above, the diffusion process was 
reduced to the solution of a single equation for the mixture fraction { and all other properties could 
then be determined by using equilibrium assumptions. For more realistic combustion problems with 
finite-rate kinetics, detailed reaction mechanisms need to be incorporated. However, it is well-known 
that using full finite-rate kinetic mechanism can cause significant numerical stiffness (due to large 
production/destruction terms) and can significantly increase the computational cost. An approach 
that reduces the computational effort is to use reduced mechanisms that result in the reduction of the 
number of equations that must be solved. Numerical computations using reduced mechanisms have 
been employed for engineering analysis (e g, Correa and Shyy, 1987; Chen and Kollman, 1990) and is 
an acceptable approach. Since our intention is to couple this mixing model to conventional prediction 
methods for engineering analysis, the demonstration of a reduced mechanism within theis mixing 
model is very important. This has been carried out for the case of a hydrogen-air nonpremixed flame 
problem. 

The hydrogen-air flame problem has been in the past both experimentally (e.g., Magre and Dibble, 
1988) and numerically (e.g., Chen and Kollmann, 1990). Also, a reduced mechanism model has been 
developed by Chen. This model has been used in the framework of the pdf method to predict the 
flow properties (e.g., Chen and Kollmann, 1990). The primary advantage of the reduced mechanism 
approach is the reduction in the number of scalars that need to be solved. For example, for the 
hydrogen-air problem, only two scalars -the mixture fraction £ and a progress variable n -need to be 
solved to obtain all other properties, as described below. Also, the reaction mechanism can be solved 
first to generate a look -up table that can then be used to interpolate for the species information 
whenever required. This clearly reduces the amount of computation required. This reduction of com- 
putational effort is important in the present case since the present model accurately computes the mix- 
ing aspect, and if the look-up table procedure can be coupled with the mixing model, then the overall 
computational cost will not be significantly more than that for the case of mixing alone. This has 
indeed been accomplished in this study. 

In the hydrogen-air nonpremixed flame problem considered here, the following ten elementary 
reactions are considered: 
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Reaction 


H + 0 2 = OH +0 
O + H 2 = OH + H 
OH + H 2 = H 2 0 + H 
OH + OH = H 2 0 + O 
H + H0 2 = OH + OH 
H + OH + M = H 2 0 + M 
H + H + M = H 2 + M 
O + 0 + M = 0 2 + M 
H + 0 +M = OH +0 
H +0 2 + M = H0 2 + M 

Here, the rate constants for the ith reaction are gi 


A i 

b. 

E, 


1.2 x !0 it 

-0.91 

69.1 

(Rl) 

1.5 xlO 7 

2.0 

31.6 

(R2) 

1.0x10* 

1.6 

13.8 

<R3) 

3.4 xlO 13 

0.0 

21.0 

(R4) 

1.5 xlO 14 

0.0 

4.2 

(R5) 

2.15X10 22 

-2.0 

0.0 

(R6) 

1.83 xlO 1 * 

-1.0 

0.0 

(R7) 

2.86 xlO 1 * 

-1.0 

0.0 

(R8) 

6.2 xlO 16 

-0.6 

0.0 

(R9) 

2.0 xlO 18 

-0.8 

0.0 

(RIO) 


by 


ki 


/t,T 4, exp 


Ej 

RTi 


(20) 


and the coefTicients 4i, and £, are obtained from Warnatz (1984) with units of cm, mol, KJ, and 
K. The third-body collision efficiencies, z,, relative to H 2 are also taken from Warnatz (1984) as z Hi 
= 1 , Zq, = 0.4, z /, 3 = 0.4, z Ui0 * 6.5, and z Ar = 0.35. The total number of scalars in this reaction 
mechanism is ten, which are the seven active chemical species (i.e., H 2 , 0 2 , H 2 0, O, H, OH, H0 2 ) plus 
temperature, pressure, and density. If we assume equal diffusivities for species and enthalpy, the mix- 
ing process can be described by a conserved scalar $ The mixture fraction is defined as the normal- 
ized mass fraction of an atomic species originating in the fuel stream. With this definition, the concen- 
trations of atomic species are linearly related to £ For example, for the H 2 -air system, three conserva- 
tion equations for H and O atoms can be derived as 


+ Hon + 2«i/ 2 + 2«ff 2 o + n HOi 


= 2 




no + n 0 a + 2 noj + n B3 o + 


2 40-Q 

Mn 3 + Wo, 


( 21 ) 

( 22 ) 


h = $ti + (1 - Qh 2 


(23) 


where n, and M, are the number of moles per unit mass and the molecular weight of the ith species, 
respectively, and <f> is the ratio between the O and N atom concentrations in air. Thus, 
<)> = 0.2 1 /0.79, and T, is the mass fraction of the ith species in the fuel. 

The partial equilibrium model was used by Janicka and Kollmann (1982) and assumes that the 
reactions R1-R5 are relatively fast compared to the three-body recombination steps R6-R10. Thus, 
the partial equilibrium model for R1-R5 provides four independent relations of the form 
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= olf n q 

(24a) 

n(pH 2 

— K-2 n OH n H 

(24b) 


— K -3 n H 7 n H 

(24c) 


— K -A n OH n OH 

(24d) 


which results in the global one-step reduced reaction mechanism 

2H 2 + 0 2 = 2 H 2 0 (25) 

The reaction rate for the global reaction can be written in terms of the rates for the elementary steps: 

U = r e + r 7 + r a + r 0 + r l0 (26) 

where r, denotes the forward reaction rate for the /th elementary step with units of moles/(mass-timc). 
With these relations, in addition to the mixture fraction £, only one more reactive scalar is required to 
describe the combustion process. In this study, following the work by Chen and Kollmann (1991), the 
total number of moles per unit mass, n, is chosen as the reactive scalar. By definition. 


n — ny 2 + ti( ^ + n y 2 o + no + n 0 y + ny + rty qj + ny 2 (27) 

The allowable domain for n at a given ( needs to be determined. The allowable maximum occurs 
when fuel and air are mixed without reaction. Then the maximum n is 


n 


max 


My 2 My 2 + 


(28) 


and the allowable minimum occurs when hydrogen and oxygen combine to form water without form- 
ing any radical. Two possible states can occur depending upon the value of the mixture fraction £ 
When £ < the minimum value is given as 


I (1-6(1+*) 

2My t M Hj + 


For £ > £* the minimum value is 


(29) 


^min 


{Yy, 

2My 3 


+ 


0-0 
+ 4>M t>) 


(30) 


Two types of hydrogen-air combustion process has been studied in this Phase I. A pure hydrogen-air 
case with a stoichiometric conditions, £* = 0.028, n H = 0.0412 and a 22% Argon+78% hydrogen -air 
case with stoichiometric conditions, = 0.1632, n* = 0.0388. The latter combustion problem was 
studied due to the availability of some experimental data (Magre and Dible, 1988) with which com- 
parison could then be made. 


A nondimensional reaction progress variable n * can be also be defined as 
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~ n 


( 31 ) 


* ftmax 

n = ^ 

^max ^min 

which has two states: n * = 0 when the mixture is completely unburnt, and n * = 1 when the mixture 
is completely burnt to form water. Then, according to the global reaction, equation (25), the rate 
equation for n * at a fixed mixture fraction £ is 

<^f>< - . ' - ['. + r, + r. + r. + r„l (32) 

ut H max ft min ' 

Either this equation or an equation for ^jj- can be used for the progress variable. In this study, the 

equation for was used. 

at 

The mechanism used to model the formation of NO in the H 2 -air combustion process is primarily 
from the thermal pathways. Thus, for now, the prompt NO formation and the formation of NO 
involving intermediate N 2 0 are neglected. This is considered acceptable for now, since in general the 
majority of NO formed is from the thermal NOx pathways. More complete NO production mechan- 
isms will be considered in the next phase. 

The thermal NO reactions are described by the well-known Zeldovich mechanism as 


*2 

+ 0 

= NO 

+ 

N 

o 2 

+ N 

= NO 

+ 

0 

N + 

OH 

= NO 

+ 

H 


(NR1) 

(NR2) 

(NR3) 


If a steady-state assumption for the N atom is assumed, and [NO]/[NO] e ,j « 1 is assumed, the Zcl- 
dovich reaction mechanism can be reduced to a single global step as 


N 2 + 0 2 = 2 NO (NR4) 

with r „ 4 = r„ u and the NO formation rate (gm/cm 3 - s) can be approximated as 

SnO a 2fc/,nrl |oJa/^o (33) 

where [C] denotes the concentration of species C (moles/cm 3 ), M NO * s the molecular weight of NO, 
andfc/ iflrl = 1.84x 10 14 exp(-38370/T) cm 3 /mole - s. 

The diffusion problem analogous to equation (1) for this hydrogen-air case is the solution of two 
equations for the mixture fraction £ and the progress variable n. They can be written as 


1341 

dt ' dx 2 


and 
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dn 

dt 


+ 


t 


( 35 ) 


= D , 


d 2 n 
dx 3 




Here, the source term in equation (35) (i.e., dn/di) can be obtained from the look-up table for the 
reduced mechanism as described above. Also, for this study, we assumed that D ( = D n = v/Sc, where 
v is determined using the Sutherland’s law. 
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4. IMPLEMENTATION OF THE MIXING AND REACTION ALGORITHMS 


The implementation of the linear-eddy model as a numerical algorithm is relatively straightfor- 
ward. However, since the description given in Section 3 covers various mechanisms that must be cou- 
pled together, the steps involved in the actual implementation are described here. Figure 5 illustrates 
the processes involved in the linear-eddy algorithm. 

The first step is to initialize the fuel flow entering the mixing zone (at x = 0) by prescribing the 
velocity, initial temperature, species distribution, Reynolds and Schmidt numbers, and so on. In a full 
implementation to the RQL combustor (to be carried out in Phase II), this initial flow field informa- 
tion will be exit plane data from the steady-state calculations in the primary fuel-rich combustion 
zone. However, if this model is used also in the primary stage (this is possible since this mixing modci 
is configuration-invariant), then the inflow conditions can be the conditions at the entrance to the 
combustor. The algorithm described in Section 3 involves six separate mechanisms: the fuel feed into 
the mixing zone, the air entrainment process throughout the mixing zone, the molecular diffusion pro- 
cess, the turbulent stirring and the finite -rate kinetics process within each cell, and the displacement 
process resulting from the effects of entrainment and thermal expansion due to combustion (discussed 
below). The implementation of these mechanisms is outlined below. 

First, molecular diffusion is carried out by time integration of the diffusion equation within the 
mixing region using a standard finite-difference approach. Typically, the time-step r for the numeri- 
cal integration is r = Cr? /D^ where n is the size of the smallest (Kolmogorov) eddy (in the imple- 
mentation, r) = Ax-nun) and C < 1 for stability. 

The other main mechanism, i.e., the triplet map event that models small-scale turbulent stirring, is 
also carried out. The triplet map does not take place continuously while the numerical integration of 
Equation (4) (or Equation (17) for the propane-air or Equations (34) and (35) for the hydrogen-air 
nonpremixed flame cases) is being carried out; rather, it occurs whenever the epoch of the inversion 
event is reached during the time integration. Typically, the epoch of the mapping event is randomly 
selected based on the overall rate, R - XL, of such events. When this time is reached, the location of 
the eddy is selected, with uniform likelihood within the cell, and a eddy size is chosen by randomly 
sampling the pdf /(/). Then the mapping process is carried cut as described earlier. This exchange is 
assumed to occur simultaneously and instantaneously so that no time elapses during this process. 
After this event, the epoch of the next event is determined based on the new value of the event rate R. 
In many cases, the time scale for the mapping is much smaller than the diffusion time step and thus, a 
series of mapping are carried out till such time as the next diffusion step is required. 

Similar to the turbulent stirring events, the entrainment of the air parcels into each cell is also 
determined based on the entrainment law and the entrainment frequency X E ; also, the epochs and 
locations of the entrainment events parallel the above-described inversion process. As the fuel and 
entrained air are being diffused and stirred within each cell of the linear-eddy, chemical kinetics will 
also be occurring that will further change the chemical composition, temperature and density in each 
cell. 

The fuel feed, air entrainment, and thermal expansion due to combustion all induce streamwisc 
fluid displacement, which represents the streamwise flow. This displacement process is implemented 
as a distinct computational step that incorporates all three contributions. After the fuel feed, air 
entrainment, and combustion during a time-step have been completed, the deviation of the cell densi- 
ties from their equilibrium values is used to determine new equilibrium values using perfect gas law 
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and a streamwise adjustment of the fluid properties in the cells carried out to adjust the properties in 
each cell. 

This completes the description of the procedure for implementing the linear-eddy algorithm to 
study unsteady mixing and combustion in a finite domain characteristic of a mixing zone. Note that 
the above discussion has identified various length and time scales that are involved in the mixing and 
combustion processes. The relationship between these scales will directly determine the effectiveness 
of the mixing strategy. 
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5. RESULTS AND DISCUSSION 

In this section, we describe the results of the various tests and simulations carried out to demon- 
strate the capability of the mixing model that was developed in this Phase I study. Since detailed data 
for the RQL combustor are not available, we decided to carry out a series of calculations designed to 
be more general than the mixing configuration representative of the RQL quench zone. This was con- 
sidered acceptable for the Phase I study since the basic technical objective was to develop a general- 
purpose mixing model that can be used in conjunction with predictive schemes. The actual imple- 
mentation of the mixing model in a practical design analysis code will be carried out in Phase II. 
Direct comparison and parametric studies of the mixing strategy for the problems of specific interest 
will also be carried out in Phase II. 

The results of the Phase I study fall into three different categories in an ascending order of com- 
plexity. The first set of calculations was carried out to demonstrate the model’s capability to handle 
multispecies mixing. Typically, the primary fuel entering the mixing zone was assumed to consist of 
five species with initially specified mass fraction distributions, and the entrained air was assumed to 
consist of two species (N 2 and 0 2 ). However, the code was written so that an arbitrary number of 
species could be simulated if required. In some of the mixing studies, the species responsible for NOx 
formation (see reactions NRI-NR4) were allowed to mix in the mixing zone, and an estimation of the 
NO production based on a global reaction mechanism was made. 

The second set of calculations was carried out for a propane-air buoyant nonpremixed flame, for 
which a simple model has been demonstrated before (Kerstcin, 1991b). In this model, all flow proper- 
ties are determined by the evaluation of a conserved scalar £ through the use of equilibrium assump- 
tions. The formation of NOx in this nonpremixed flame problem was also estimated using a global 
mechanism as described in Section 3.3.2. 

Finally, the third set of calculations was closer to the more representative combustion approach 
used in many calculations for practical applications (e.g., Chen and Kollmann, 1990, 1991). For 
demonstration, we chose to study the H 2 -air combustion system. To reduce the complexity of the 
finite-rate kinetics and to reduce the computational cost, a reduced-mechanism approach was 
employed. A look-up table for the H 2 -air combustion was created in terms of two scalars: a con- 
served scalar £ and a progress variable n\ An option was included to generate a look-up table in terms 
of three scalars (i.e., £ n, and h), where h is the enthalpy of the mixture. This option is typically used 
when the effect of radiation heat loss needs to be incorporated. However, for the Phase I study, it was 
decided to use only the two-scalar look-up table. This look-up table was incorporated within the 
framework of the mixing model, and a series of computations was carried out for the H 2 -air 
nonpremixed flame. 

Since the computations are carried out for a long time, the spatial distribution will change from 
instant to instant. As the time increases, the flow eventually reaches a statistically steady state; the 
instantaneous realizations can be Favre-averaged to obtain the steady picture of the flow field. To 
obtain the steady-state field, the total number of grid points (1250 in the present case) is divided into 
blocks each containing 50 cells; each block is then further divided into bins containing 10 cells each. 
Favre-averaging of the flow variables is then carried out over each bin (10 cells) to obtain the aver- 
aged property. In general, the Favre-averaging process results in 
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( 36 ) 
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where <j> denotes any flow property of interest, V is the volume of each cell, and the summation is for 
each bin containing 10 cells. Most of the results presented in this section are in terms of these Favre- 
averaged properties, since these are typically the properties required for design and prediction. How- 
ever, note that since the integral length scale, L(x), and Ax are both increasing in the streamwisc 
direction, the total volume in each block is also increasing in the streamwise direction. Thus, when the 
Favre-averaged properties are plotted as a function of the blocks, in the physical space (i.e., x /d 0 ) the 
scale is nonlinear. This is reflected in most of the figures as noted below. 

All computations were carried out on an IBM PC compatible (486), since access to a supercom- 
puter was not provided in the Phase I study. The resource and time constraints of the Phase I study 
limited the number of the simulations that could be carried out. 

5.1 Multispecies Mixing and NOx Formation 

The first series of test cases involved modeling the effect of mixing between multispecies. For 
these tests, we assumed that the primary fuel stream consists of five species (0 2 , N 2 , O, N, and NO) 
with specified mass fractions (F.) and with an initial temperature of 2100 K. The entrained air was 
modeled by two species (0 2 and N 2 ) at 300 K. The reference primary fuel jet speed was 12 m/s and 
the secondary air was initially at stagnant conditions. Tests to study the effects of varying the fuel jet 
speed, the initial mass fraction distribution, and the entrainment conditions were carried out. Only 
representative results are shown here. Table 1 shows the various test cases discussed in this section. 

Figure 6 shows a test case (M-l) with turbulent mixing and air entrainment occurring for the full 
extent of the jet. A spreading parameter of c = 0.05 was used in equation (6), which essentially modi- 
fies the air entrainment process. Also, since now L(x) does not increase rapidly in the downstream 
direction, for the resolution used (1250 cells) the streamwise extent for these simulations was limited 
to xfd 0 < 50. Figures 6a through 6d show, respectively, the species mass fractions Y N3 , Y NO> and Y 0 
and the mixture density p/p„. Since this simulation is a pure mixing case, there is no reaction occur- 
ring and the flow field shown here is representative of the single realization of the axial distribution of 
the flow properties. Figure 6a shows that near the exit of the primary fuel jet there is hardly any 
nitrogen present (due to the initial condition, see Table I), but as the air is entrained, the nitrogen con- 
tent increases downstream. The axial variation of all the species looks the same since all species diffu- 
sion coefficients were assumed to be the same, and the effect of turbulent mixing on all the species 
was also the same. As the jet spreads due to air entrainment, the density (Figure 6d) increases in the 
far field, since there the mixture is mostly air. 

Figure 7 shows another test case which was similar to case M-l shown in Figure 6 except that for 
x/d 0 < 1 3 (or the first ten blocks), no stirring (turbulent mixing) or air entrainment was allowed to 
occur. Constraints such as this are one way to mimic different mixing strategies. For this condition, 
the constraints imply that for x/d 0 < 13, the mixture from the inlet undergoes simple diffusion without 
mixing with the secondary air. This is clearly reflected in the instantaneous species, density, and tem- 
perature distribution in the jet as shown in Figure 7. Also, in this case, NO production based on the 
global hydrogen-air jet flame model (Chen and Kollmann, 1991) given earlier was used to estimate the 
NO production (test case M-3). The results of three cases (M-l through M-3) are summarized in 
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Figures 8a through 8c, which show the Favre-averaged flow properties. Note that in these plots the 
x-axis is in terms of the blocks (1 to 25), which is related to the streamwise distance ( x/d 0 ) non- 
linearly as shown in Figure 8a. As noted before, as the block location increases downstream, the phy- 
sical spacing increases due to the increase in the cell size Ax. For case M-3 (with NO production), it 
can be seen that the peak NO production is close to the location where air was allowed to be 
entrained ( x/d 0 = 13, Figure 8b). This is a consequence of suddenly allowing the air to be entrained 
in a large parcel (since the volume of the block is large), which then mixes with the hot primary 
stream causing rapid production of NO. Since NO is being produced at the expense of O-atoms, the 
mass fraction for oxygen atoms decreases rapidly (Figure 8c) for test case M-3. As the O atoms 
become negligible, the NO production also decreases very rapidly. Note that, for this test case, there 
are sufficient N 2 molecules available for the NO production mechanism. Thus, for this case, the 
important species is the O-atom, which was only available from the initial distribution from the fuel 
jet. 

Parametric studies such as this with various entrainment conditions may be able to identify the 
important mixing strategy that will minimize NO production in the RQL combustor. 

Figure 9 shows another test case (M-4) for which the initial species distribution in the fuel jet was 
adjusted so that the fuel stream entering the mixing zone contains no N 2 molecules and no NO 
molecules (see Table I). Also, in this case, the free jet spreading rate c =0.14 was used, and mixing 
and air entrainment were allowed to occur the entire extent of the mixing zone. In these calculations, 
the streamwise length of the mixing region exceeds x/d 0 = 200. Figure 9a shows the spatial distribu- 
tion of for this test case. As the jet spreads, around x/d 0 = 25, the entrainment of air (which 
contains a significant amount of N 2 ) causes the increase of N 2 content in the mixing region. Mixing 
with the cold air stream causes the temperature of the fuel jet to decrease very rapidly. A conse- 
quence of this effect is that NO production is limited to a small region near the zone where the 
entrained air is mixing with the fuel stream. This is seen in Figure 9b. Figure 9c shows the Favre- 
averaged distribution of the N 2 , 0 2 and NO molecules as a function of the blocks. Again, peak pro- 
duction of NO is observed in the region where the hot fuel stream is mixing with the cold entrained 
air at around x /d 0 = 25. 

Finally, for reference, we carried out a simulation using representative conditions as reported in 
Howe et al., (1991) with the conditions at the exit of the primary stage of the RQL combustor entering 
the quench zone (see Table I). Figures 10a through 10c show, respectively, an instantaneous realiza- 
tion of the temperature, the oxygen mass fraction Y ^ and the water mass fraction Y u 2 o variation in 
the streamwise direction for this test case. No combustion was allowed to occur, and this simulation 
was carried out just to see the extent to which the streamwise penetration of the fuel stream will occur 
in the absence of reactions. 

5.2 Propane-Air Diffusion Flame 

A series of test cases was studied similar to the cases studied earlier by Kerstein (1991a). This is a 
model problem of a buoyant propane-air nonpremixed flame. Both the primary fuel and the secon- 
dary air were assumed to be at the same temperature (300 K) and the flame temperature was 2280 K 
(see Table II). The effect of increasing the fuel jet velocity was studied. Figure 1 1 shows the typical 
instantaneous realization of the axial variation of the various flow properties for a specified fuel jet 
speed (test case P-1). As the fuel enters the mixing region, whereever the mixture fraction becomes 
stoichiometric, the fuel is converted to product and the temperature rises to the adiabatic flame 
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temperature. This is reflected in Figures 1 la and 1 lb. Further downstream only air is present, which 
is at a much higher density. 

Figure 12 shows the instantaneous flow properties for another test case (P-4) which was for a much 
higher jet speed of 72 m/s. In this case the fuel penetrates further downstream in the streamwise 
direction (Figure 1 2a), and more regions with high temperature in the mixing region are seen (Figure 
12b). 

Figure 1 3a shows the time-averaged density field in the jet for various jet velocities. As the fuel 
jet speed increases, the increase in density occurs further downstream. This is similar to the results 
obtained earlier by Kerstein (1991a). Figure 13b shows the effect of carrying out the computations 
twice as long. Here, / x »s the time at which the simulation was stopped. Typically, the simulations 
should be carried out long enough so that the flow properties reach a statistically steady state. Figure 
1 3b shows that, for this simulation, the far field has still not reached a statistically steady state; how- 
ever, the difference between the two simulations is not very significant. Resource and time con- 
straints of the Phase I study made it difficult to carry out the simulations long enough to reach a sta- 
tistical steady state everywhere. The intent of this Phase I study was only to demonstrate the potential 
capability of the mixing model. 

Validation of the mixing model was carried out earlier by Kerstein (1991a) where the simulations 
were performed for the propane-air case in more detail. Since that study is quite relevant for the 
present work, the important conclusions are briefly noted here by using some of those results. Figure 
13c shows the time-averaged density field as a function of streamwise distance (which is non- 
dimensionalized by the jet momentum diameter, d e ). The linear-eddy model predictions are com- 
pared to experimental results of Becker and Yamazaki (1978) and to the predictions of an "instant 1 ' 
mixing model. Clearly, the linear-eddy model is capable of providing improved predictions when 
compared to the instant mixing model. 

Since NOx production was allowed to occur based on the global mechanism as given in equation 
(19), it was possible to estimate the emission index for NO as a function of time of evolution for the 
entire mixing region. Figures 14a through 14c show the NO emission index as a function of nondi- 
mcnsional time for various fuel jet speeds (NOTE: the emission index is given here in kg /kg of fuel 
and to determine the standard value in gm /kg of fuel the value igiven here must be multiplied by 
1000). Here, the time t = t /t 0 , where t 0 = d e /u itt . Thus, for increasing jet speeds, t 0 becomes smaller. 
We must point out that the NO emission estimate shown in Figure 14 was obtained only at a few times 
(marked by solid circles), and the large fluctuation in time seen in these figures may not be meaningful 
since the in-between time values were not recorded. Also, the simulations were not carried out long 
enough to reach a statistically steady state due to resource (computer) and time constraints. 

Figures 1 5a through 1 5c show the axial variation of the NO production rate for various jet speeds. 
These plots indicate that as the jet speed increases the peak production of NO moves further down- 
stream. Figure 1 5c also shows the effect of carrying out the simulation for twice as long. The flow 
field has clearly reached a statistically steady state at nearly the entire jet except for the region of 
peak NO production. 

5.3 Hydrogen-Air Diffusion Flame 

For this test case, we chose to specify the initial values of the conserved scalar £ and the progress 
variable n in the fuel jet to represent the H 2 fuel. Thus, the fuel jet enters the mixing domain with 
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£ = 1 and n = 0.496 (for the 22%Ar+78%hydrogen case, n = 0.0965) and with a temperature of 300 
K. The stoichiometric conditions for the H 3 -air combustion are given by £,« = 0.028 and n * = 0.0412, 
and the flame temperature is 2193 K. For the 22%Ar+78%H 3 -air combustion case, the stoichiometric 
conditions are given by £* = 0.1636 and n* = 0.0388 and the flame temperature is 2303 K. The air 
stream entrained into the mixing zone was assumed to be at a state given by £ = 0 and n = 0.03466, 
and a temperature of 300 K. The effects of changing the fuel jet diameter, the fuel jet velocity, the *» 
spreading (and hence the entrainment) rate, and the location of entrainment were studied; only 
representative results are shown here (see Table III). 

Figures 16a through 16d show, respectively, an instantaneous realization of the axial distribution 
of the scalar £, the progress variable n, the NO moles njvofboth n and « wo are in moles/unit mass) and 
the mixture density p (in kg/m 3 ) for a representative test case. For this simulation, a 0.52-cm diame- 
ter fuel jet enters the mixing region at a jet speed of 1 50 m/s and was modeled with the entrainment 
of air occurring over the entire axial extent with the free-jet spreading rate noted earlier. Figure 16a 
shows that the conserved scalar £ starts with a value of 1 from the fuel jet (at x /d 0 = 0) but rapidly 
drops to low values as the air is entrained and combustion occurs. The progress variable n also shows 
a similar variation. Since both the mixture fraction and the progress variable were diffused and mixed 
identically, the variation in space at any given time is the same. The instantaneous variation of the 
NO moles, n NO , shows that most of the local NO is produced around x /d 0 = 65. Finally, Figure 16d 
shows the density field at a given instant, indicating an increase in the density near the outflow. This 
increase is mainly due to the fact that most of the flow far downstream primarily contains air, which 
has a much higher density than H 3 . 

The axial variation of the species in the H 3 -air reaction mechanism can also be determined for 
each of these test cases. Figures 17a through 17c show the axial variation (in moles/unit mass) of the 
species H 3 0, H 3 , 0 3 , H, O, OH, and H0 3 for a simulation which was similar to that shown in Figure 
16 except that the fuel jet speed was 75 m/s. This configuration is similar to some experiments car- 
ried out in the past at Sandia Laboratories. To understand the mechanism of NO production and the 
variation of the different species in the mixing region, each of the parts of Figure 17 shows multiple 
plots for the various flow properties. For example. Figure 17a shows the temperature and the concen- 
trations of O-atoms and NO molecules in the mixing region. Clearly, peak NO production occurs at 
the same location where temperature is also at a maximum. Since NO production requires the pres- 
ence of O-atoms by virtue of the production mechanism, most of the NO is concentrated at the loca- 
tions where O-atoms are present. As O-atoms are depleted, the NO production also decreases rapidly. 
Figure 17b shows how the other flow properties, for example, the concentration of H 3 0, H 3 , 0 3 , and 
OH, vary with axial distance. As combustion occurs, the H 3 concentration is depleted while the 
major product (water) begins to form and the temperature rises. The oxygen concentration increases 
in the far field due to the entrained air. The peak production of OH occurs close to the location of 
peak temperature. Thus, it appears all the radicals are located close to the peak temperature region. 
Finally, Figure 17c shows the concentrations of the rest of the species, H, OH, and H0 3 for this simu- 
lation. As the H 3 in the fuel jet breaks down H-atoms are formed much earlier; some of the H-atoms 
are combined with the available O-atoms to form OH near the peak temperature location. However, 
some of the H atoms are transported further downstream due to turbulent mixing to form H0 2 as can 
be seen in this figure. 

Figures 1 8a and 1 8b show another simulation similar to that shown in Figure 1 7 except that, in this 
case, the fuel jet speed is doubled to 150 m/s. All conclusions discussed above for the lower jet speed 
case (Figure 17) hold again for this simulation. The only major difference is that the peak locations 
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appear to move slightly further downstream. To avoid causing confusion (due to the presence of mul- 
tiple curves), the scales corresponding to the various flow properties are not shown in these two fig- 
ures. 

Figures 19a and 19b show the axial variation of the NO emission index for these two calculations. 
The emission index starts to increase at around x/d 0 = 40 for both cases and in the far field levels off. 
A similar trend has been observed in earlier H 2 -air experiments. Note that with an increase in jet 
speed the emission index decreases. This also agrees with earlier experimental data (Chen and 
Driscoll, 1991), although a direct comparison is not possible since the conditions used in the present 
test case and those used in experiments are not the same. Furthermore, the present simulation has not 
been carried out for a sufficient long time to reach a statistically steady state. This can be seen 
because the length of the flame can be estimated from Figure 17a, by locating the downstream dis- 
tance at which the mixture fraction becomes stoichiometric. The present simulation indicates a flame 
length of around x /d 0 — 75, which is less than the experimentally observed value. 

Figures 20a and 20b show respectively, the scatter plot representation of T—( variation with time 
at a given axial location of x/d 0 = 30. The general trend is very similar to earlier experiments as will 
be described further below. 

Finally, a series of calculations was carried out using the 22%Ar-78%hydrogen flame combustion 
so that comparisons with earlier experiments (Magre and Dibble, 1988) and pdf computations (Chen 
and Kollman, 1990) could be carried out. Figures 25a through 25c show, respectively, the present 
computations, the experimental data of Magre and Dibble (1988), and the pdf computations of Chen 
and Kollmann (1990) in terms of the T-£ scatter at x/d 0 = 30 for fuel jet speed of 75 m / s. Figures 
26 and 27 show, respectively, these three plots for the jet speeds of 150 and 225 m/s. Clearly, the 
present linear-eddy computations are capable of reproducing quite accurately the experimental and 
earlier pdf computed predictions. The increased scatter in the experimental data is due to the noise in 
the measurement, but the agreement between the present computations and these earlier works is 
quite remarkable. 

Figure 24a shows a three-dimensional surface plot of the NO moles/unit mass n NO as a function of 
the mixture fraction £ and progress variable n for test case H-4. This figure clearly shows the allow- 
able domain in the «-£ space where NO production occurs. It can be seen that peak production of 
NO occurs near stoichiometric conditions. Figure 24b shows another three-dimensional NO surface 
in the 7"-£ plane. NO is primarily produced at locations near stoichiometric conditions and at high 
temperatures. Finally, Figure 24c shows the variation of the NO Emission Index as a function of fuel 
jet speed. Also shown in this figure is the recent pdf prediction of Chen and Kollmann (1991). The 
linear-eddy prediction also shows the same trend, i.e., the NO Emission Index decreases with increas- 
ing jet velocity. However, the current predicted NO Emission Index does not agree with the pdf 
results for several reasons. The present calculation has not been carried out long enough to reach a 
statistically stationary state and, therefore, the present results do not indicate the final state in the 
mixture. Also, lack of time and resources also made it difficult to evaluate if the resolution used in 
the present calculation was sufficient. The pdf computations employed a surrounding coflowing air at 
a velocity of around 9.2 m/s (the linear-eddy calculations assumes that there is no coflowing air and 
the entrained air is initially at stagnant conditions). Finally, the specific geometry used in the pdf cal- 
culations was different (coaxial tubes) from the geometry modeled in the present study. All of these 
factors could contribute to the observed differences between the current computations and the past 
studies. However, the ability of the present mixing model to predict experimentally observed trends 
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(i.e., the decrease in emission index with jet speed) is very encouraging. 

Finally, to demonstrate that the present model can predict experimentally observed values reason- 
ably well. Figure 25 shows an earlier computation by Kerstein (1991a) for the H 2 -air nonpremixed 
flame that was compared to some earlier data from Sandia. The reaction mechanism for this study 
was not the reduced-mechanism used here; an equilibrium assumption was used to relate all properties 
to the mixture fraction, as was done in the propane-air nonpremixed flame study described in this 
report. Clearly the capability of this mixing model to study realistic problems and its superiority over 
the instant mixing model have been established by the present study and by the earlier studies by 
Kerstein. 
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6. CONCLUSIONS AND RECOMMENDATIONS 


In this report, the development of a new mixing model is described that takes into account the fun- 
damental features of turbulent mixing. Both turbulent mixing and molecular diffusion in the small 
scales have been identified in experiments to be important in flows with chemical species, fey specifi- 
cally separating these two effects and modeling them in the small scales, the mixing model developed 
here has a capability that all previous mixing models have lacked. This mixing model was used to 
study multispecies mixing, a propane-air nonpremixed flame, and a hydrogen-air nonpremixed flame. 
Since these three test cases required different model formulations, the capability of the mixing model 
to handle all three different test cases clearly show its potential for wide application. The production 
of NO due to a thermal mechanism (the Zeldovich mechanism) was also estimated in all three cases, 
and the effects of jet diameter, jet speed, and modification of the turbulent mixing and air entrain- 
ment process were studied. Due to resource limitations, a parametric study has not yet been carried 
out; however, the series of calculations shown in this report indicates that a parametric study can 
easily be accomplished. 

Comparison of the present results with experimental data and the earlier comparisons by Kerstein 
(1991a) clearly demonstrate the strengths of the linear-eddy model. The present study has shown that 
finite-rate kinetics in the form of a reduced-mechanism can easily be incorporated within the frame- 
work of the mixing model. The work carried out so far has provided enough understanding of the 
behavior of the mixing model and has laid the basic framework for application to practical design 
problems, which will be studied in the next phase. 

In the next phase (a proposal will be submitted), we propose to further improve the model and 
implement it in a manner that would provide a numerical capability to address practical issues such as 
the mixing process in the quench zone of the RQL combustor. Issues regarding the grid resolution, 
the time of simulation, and the mixing geometry will be addressed first. Then, a series of parametric 
studies will be carried out with the present code to determine the importance of the various time and 
length scales that were identified in the Phase I study. 

In parallel to this parametric study, a major research effort will be undertaken in which the present 
mixing model will be coupled to a steady-state prediction code. The steady-state code will not be 
developed. Two possible options are currently viable. The first approach is to couple the mixing 
model to a standard finite-difference steady-state code using the k - t turbulence model. The second 
. oproach is to couple the present mixing model to a prediction code based on the pdf method for the 
joint-scalar pdf or the joint-velocity-scale pdf. The linear-eddy model has many elements that are also 
present in the Monte Carlo approach used in pdf methods. Therefore, it is possible to implement this 
model within the pdf method to account for the unsteady mixing effects. The second approach (cou- 
pling the mixing model to the pdf method) may be, in the long run, more productive, since the pdf 
method is sup-*, : r to the standard finite-difference method. The basic weakness of the current pdf 
codes is in the mixing model (typically, mixing models similar to Curl’s mixing model is used) and, 
therefore, the strength of the present linear-eddy model could be used to enhance the pdf method’s 
capability. Of course, to couple the present mixing model to the pdf method will require access to the 
code. Pdf codes are used in General Electric and General Motors research groups (private communi- 
cations); however, these codes are not in the public domain. We will discuss the issues involved in 
coupling the mixing model to steady-state codes in more detail in the Phase II proposal. 
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Implementation of the present model within the framework of a standard finite-difTerence code or 
a pdf code will require significant modifications to the present model (particularly the application- 
specific elements described in Section 3.2). Fcr example, if this mixing model is used in a standard 
finite-difference code, then instead of using the discretization shown in Figure 2 for the entire mixing 


region, the mixing model will be implemented within each grid cell as shown in Figure 26. Computa- 
tions of turbulent mixing, molecular diffusion and chemical heat release can then be carried out 
within each cell as described in Sections 3.1 and 3.3. Large-scale transport due to the mean flow and 
turbulent fluctuations primarily due to the volume fluxes from the neighboring cells, will result in 
changes to the flow properties in each cell. The exchange of species between neighboring cells will 
have to be taken into account. All these processes will result in a new state of the species, tempera- 
ture, and density field within each grid cell. By Favre-avcraging these properties within the linear 
eddies in each cell, the new mean properties can be determined and can then be used to correct the 


mean mass, momentum, and energy transport being computed in the steady-state algorithm. A tech- 
nique for modeling the exchange between neighboring cells and a method for coupling the effects of 
small scale mixing and combustion to the Navier-Stokes equations (mass, momentum, and energy 
equations) have recently been developed (Menon, 1991; Menon et al., 1991; McMurtry et al., 1992) 
within the framework of a subgrid model approach. Therefore, the specific method for coupling the 
present model to standard codes has already been developed. The coupling issues will be will be dis- 
cussed in more detail in the Phase II proposal. 


Finally, the physical problem discussed in this report is the mixing process and NO formation in 
the quench zone of the RQL combustor. However, from the formulation described in this report, it 
should be clear that the present model is not restricted to that mixing geometry and could be used lor 
other flows. If the quench zone of the RQL combustor is modeled, then to study the mixing problem, 
the inflow to the mixing model (i.e., the exit conditions from the primary stage) must be provided from 
some steady-state prediction. However, if the original steady-state prediction of the state of the flow 
in the primary stage is in error (a serious possibility), then even if the present model is accurate, its 
predictions will be contaminated by the errors in the original input. A more practical approach would 
be to use the present model within the framework of a steady-state code to first predict the flow pro- 
perties in the primary stage before addressing the quench zone. All these issues will be discussed in 
more detail in the Phase II proposal. 
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Table I. Multispecies Mixing Test Cases 



Condition B: = 0.0, = 0.6, Y N = 0.1, Y Q = 03, Y NQ = 0.0 

Condition C: Y^ = 0.0, Y^ = 0.6925, Yj^ = 0.1354, Y co = 0.1578, Y^ = 0.0143 



Table II. Propane- Air Nonpremixed Flame Test Cases 



Fuel Jet Conditions 

Test 

Case 

Vel 

(m/s) 

Temp 

r/r. 

Density 

(p/p») 

Dia 

(m) 

P-1 

12 

7.61 

1.526 

0.0039 

D 

24 

7.61 

1.526 

0.0039 

D 

48 

7.61 

1.526 

0.0039 

S 

72 

7.61 

1.526 

0.0039 


Table III. HydrogervAir Nonpremixed Flame Teet Caee 


■ 

H 

Fuel Jet 

Conditic 

ns 

Stoichiometry Conditions 

Comments 



Temp 

•K 

Dia. 

(m) 

* * O 

Initial Cond 


*• 

H 


H-l 

75 

300 

0.005 

1 -i 

r\ 

0.0412 

0.028 

2193 

L * 1 + 0.14x, H 2 *air combustion; 
reduced mechanism 

H-2 



0.005 

\ * 1 

0.0412 

0.028 

2193 

L ■ 1 + 0.14x, H 2 'air combustion; 
reduced mechanism 

H-3 

75 

300 

H 


0.0388 

0.1632 

2303 

L * 1 + 0.14x, 22% At + 78% H 2 ~air combustion; 
j reduced mechanism 

H-4 

150 

300 

0.0052 

VA 

•n -o aus 

0.0388 

0.1632 

2303 

L * 1 + 0.14x, 22% Ar + 78% H 2 *air combustion; 
reduced mechanism 

H-5 

225 

300 

0.0052 

T\ :(J.oUS 

0.0388 

H 

2303 

L * 1 + 0.14x, 22% At + 78% H 2 -air combustion; 
reduced mechanism 
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Secondary Air 



Figure 1. Idealized RQL combustor 



Figure 2. Mixing region downstream of a fuel jet 


BLOCK INVERSION 



TRIPLET MAP 


Figure 3. Alternative rearrangement processes representing the effect of individual turbulent 
eddies. In each panel, the double-arrowed line represents the spatial domain. Tick marks 
demarcate a segment selected for rearrangement. For the purpose of this illustration, the 
concentration field prior to rearrangement is taken to be linear in the spatial coordinate. 
Block inversion (a), employed in previous applications, introduces two discontinuities in the 
concentration field, indicated by the vertical dashed lines. The triplet map (b), adopted in the 
present formulation, introduces no discontinuities, though it introduces discontinuous 
derivatives. It triples the number of level crossings of any concentration value within the 
segment (e.g., crossings of the value represented by the horizontal dashed line), corresponding 
to a tripling of area of any material surface within the segment. (From Kerstein, 1991a) 
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c(x) 


X 


t 



Figure 4. Schematic illustration of the effect of a single, clockwise eddy on a two-dimensional 
scalar field that initially has a uniform concentration gradient. Upper panel: initial 

concentration isopleths (vertical lines) and concentration profile c(x) (heavy line) parallel to 
the initial concentration gradient. Lower panel: concentration isopleths and concentration 
profile at a later time. (From Kerstein, 1991a) 



Figure 5. Block schematic of the linear-eddy algorithm 
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Figure 6. Axial variation of the species and flow properties during multispecies mixing. Test 
condition M-l (see Table I). 
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8 16 


O 



a. Nitrogen mass fraction (K/^) 

b. Nitric oxide mass fraction (IVo) 

c. Oxygen atom mass fraction ( Y Q ) 

d. Density (p/p a ) 

e. Temperature ( T /T a ) 


Figure 7. Axial variation of the species and flow properties during multispecies mixing. Test 
condition M-3 (see Table I). 
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Figure 8. Axial variation of Favre-averaged flow properties for test cases M-l, M-2 and M-3 









a. Nitrogen mass fraction (K Wa ) 


b. Nitric oxide mass fraction (K no) 



0 tTs 2L75 5T0 HOO 2?8.0 

*/d 0 

c. Axial variation of Favre-averaged flow properties. 


Figure 9. Axial variation of the species and flow properties during multispecies mixing. Test 
condition M-4 (see Table I). 
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b. Water mass fraction (Y u 2 o) 


Figure 10. Axial variation of the species and flow properties during multispecies mixing. Test 
condition M-5 (see Table I). 








x /d 0 

a. Density variation with fuel jet speed 



x /d 0 


b. Density variation with time of simulation 



x/dg 


c. Comparison of computed and experimental axial density variation 
(Kerstein, 1991a) 

Figure 13. Axial variation of time-averaged non dimensional density (p/p a ) during propane-air 
combustion for various test cases. 


44 






NO EMISSION INDEX * no EMISSION INDEX 



time 0/t o ) ^ me 

a. NO Emission Index for test case P-1 (Uj e t =12 m js) b. NO Emission Index for test case P-2 (u jri = 24 m j s ) 



c. NO Emission Index for test case P-3 (u jet = 48 m j s ) d. NO Emission Index for test case P-4 (u jet = 72 m js) 


Figure 14. Evolution of the NO Emission Index (in kg /kg of fuel) as a function of time for 
various fuel jet speeds. Propane-air combustion for test cases P-1 through P-4. 

Solid circles denote locations where NO Emission Index was recorded. 


45 





NO EMISSION INDEX 


0 


7.25 


0.0 


7.25 2 1.75 x /do 5 1.0 110.0 

b. Test case P-2 (u }et = 24 m /s) 


x/d 0 51.0 nu.u 228.U 

a. Test case P-1 (it )et = 1 2 m /s) 



c. Test case P-3 ( u )ct = 48 m /s) 


Figure 15. Axial variation of nitric oxide production rate for various fuel jet speeds. Propane- 
air combustion cases P-1 to P-3. 












o.o 
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_ \ _ 

»• Nitric Oxide Oxygen atom (ii 0 ) in | kmolcs /kg] and Temperature (T) [K) 



0.0 7.25 21.75 51.0 110.0 228.0 

x/d Q 

b. Water ( n Ba0 ), Hydrogen Hydroxyl (« 0 i i)» Oxygen (hq,) in [kmoles /kg] and 

Temperature (7) in |K] . 



0 0 5.0 10.0 15.0 20,0 25 .0 

.6 725 21.75 5L0 110.0 "228.0 x /d 0 

c. Hydrogen atom Hydroxyl (n OH ), and Hydrogen dioxide (n a02 ) in [kmoles /kg\ 


Figure 17. Axial variation of Favre-ave raged species during Hydrogen-air combustion. Test 
case H-2 (see Table III). 
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Blocks 

x/d 0 


a. Mixture fraction (& progress variable (w), Water (n Hj0 ), Nitric Oxide (n jV o) in [kmoles /kg]. 
Temperature (7) in |K| and density (p) in \kg /m 3 1 



b. Hydrogen («a 2 ), Hydroxyl (h 0 h)i Oxygen («pj)> Hydrogen atom (n#), Oxygen atom 
and Hydrogen dioxide (n UOj ) in \kmoles /kg] 

Figure 18. Axial variation of Favre-averaged flow properties during Hydrogen-air combustion. 
Test case H-3(see Table III). 
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Figure 19. Axial variation of the NO Emission Index (in kg /kg of fuel) for different fuel jet 
speeds. 
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a. Test case H-3 (see Table III) 



b. Experimental data (Magre and Dibble, 1988) 



c. Pdf calculation (Chen and Kollmann, 1990) 


Figure 21. Comparison of computed and experimental scatter plot of temperature as a function 
of mixture fraction. Hydrogen-air combustion (with 22%Argon + 78%Hydrogen in the fuel 
jet) (itj't = 75 m/s) 
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T[K 1 



b. Experimental data (Magre and Dibble, 1988) 



c. Pdf calculation (Chen and Kollmann, 1990) 


Figure 22. Comparison of computed and experimental scatter plot of temperature as a function 
of mixture fraction. Hydrogen-air combustion (with 22%Argon + 78%Hydrogen in the fuel 
jet) (m jn = 150 m/s) 
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b. Experimental data (Magre and Dibble, 1988) 



c. Pdf calculation (Chen and Kollmann, 1990) 


Figure 23. Comparison of computed and experimental scatter plot of temperature as a function 
of mixture fraction. Hydrogen-air combustion (with 22% Argon + 78%Hydrogen in the fuel 
Jet) Ui,et = 225 m /. s ) 
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c. NO Emission Index (gm / kg of fuel) 
as a function ot jet speed 


Figure 24* The NO Emission dependence on temperature, mixture fraction, progress variable 
and fuel jet speed. 22%Argon+78%Hydrogen combustion case. 





x/dg 


Figure 25. Axial variation of time-averaged density. 22%Argon+78%Hydrogen-air combustion 
case modeled using equilibrium assumptions (Kerstein, 1991a)* 




Figure 26. Schematic of a spatially developing shear layer and the implementation of the 
linear-eddy mixing model within the large-scale grid. 
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